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Doubly virtual Compton scattering and the beam normal spin asymmetry 


Mikhail Gorchtein 1 

1 California Institute of Technology, Pasadena, CA 91125, I/5b0 

We construct an invariant basis for Compton scattering with two virtual photons (VVCS). The 
basis tensors are chosen to be gauge invariant and orthogonal to each other. The properties of the 
corresponding 18 invariant amplitudes are studied in detail. We consider the special case of elastic 
VVCS with the virtualities of the initial and final photons equal. The invariant basis for VVCS 
in this orthogonal form does not exist in the literature. We furthermore use this VVCS tensor for 
a calculation of the beam normal spin asymmetry in the forward kinematics. For this, we relate 
the invariant amplitudes to the helicity amplitudes of the VVCS reaction. The imaginary parts 
of the latter are related to the inclusive cross section by means of the optical theorem. We use 
the phenomenological value of the transverse cross section or ~ 0.1 mbarn and the Callan-Gross 
relation which relates the longitudinal cross section ol to the transverse one. The result of the 
calculation agrees with an existing calculation, and predicts negative values of the asymmetry B n 
of order of 5 ppm in the energy range from 3 to 45 GeV and for very forward angles. 

PACS numbers: 12.40.Nn, 13.40.Gp, 13.60.Fz, 13.60.Hb, 14.20.Dh 


I. INTRODUCTION 

Over recent years, interest in Compton scattering with 
two virtual (typically, space-like) photons has arisen and 
has been constantly growing. This interest is based on 
significant progress in precision electron-proton scater- 
ing experimental technique. Polarization transfer data 
for elastic electron-proton scattering fl] showed that the 
extracted values of the electromagnetic form factors ob¬ 
tained with this technique differ significantly from the 
Rosenbluth separation data Q ■ It has been pointed out 
that these two data sets may be reconciled by the proper 
inclusion of higher order contributions in the electromag¬ 
netic coupling constant a emi namely the two photon ex¬ 
change contribution |3|]. Several model calculations using 
the elastic form factors only 0 , the A(1232) resonance 
@, and the partonic picture ]b], show that the contri¬ 
bution of the two photon exchange exhibits the correct 
trend which may be enough to explain the discrepancy of 
the experimental results. On the other hand, a different 
class of observables has become experimentally accessi¬ 
ble in the framework of parity-violating electron-proton 
scattering. These new observables involve the spin orien¬ 
tation normal to the reaction plane and have been shown 
to be directly related to the imaginary part of the elastic 
electron-proton scattering amplitude. Since one photon 
exchange does not lead to an imaginary part, the lead¬ 
ing order contribution comes from the exchange of at 
least two photons [j]. Since the exchange of three or 
more photons is suppressed by powers of the small cou¬ 
pling constant a e m , a measurement of such an observable 
probably provides the cleanest way to access Compton 
scattering with two virtual photons. At present, several 
experimental data points exist in different kinematical 
regimes |§. On the theoretical side, significant progress 


has been made towards quantitative understanding of the 
experimental data Q. 

These experimental and theoretical studies require a 
description of Compton scattering with two virtual pho¬ 
tons. The general form of the explicitly gauge invariant 
basis for VVCS was developed in the 1970’s by Tarrach 
0 - The invariant amplitudes of Tarrach have the ad¬ 
vantage that they don’t have kinematical singularities 
nor constraints, and thus can be used for a dispersion 
calculation. However, they were never used for a realis¬ 
tic calculation. One of the reasons is that, at that time, 
the VVCS reaction was of pure theoretical interest (apart 
from its forward limit which is used in DIS and described 
by the nucleon structure functions rather than the in¬ 
variant amplitudes of Tarrach). The other reason is the 
technical difficulty in using this basis in a real calcula¬ 
tion, which is due to the non-orthogonality of the basis 
tensors. Another approach was used by Prange for real 
Compton scattering 11] and then generalized to Comp¬ 
ton scattering with a virtual initial photon by Berg and 
Lindner (l2H |. This method consists in defining the in¬ 
variant Compton basis in terms of gauge invariant basis 
tensors which are orthogonal to each other. In this work, 
we generalize this approach to the case when the outgo¬ 
ing photon is also virtual. In section [H] we introduce the 
VVCS tensor in such a form, and study the properties 
of the invariant amplitudes. We also consider the special 
case of elastic VVCS with initial and outgoing photons of 
equal virtuality. In the literature, the invariant Comp¬ 
ton tensor for this reaction exists only for the forward 
case. We obtain here the result for the general case. In 
section UTT1 we use the introduced from of the Compton 
tensor for a practical calculation of the beam normal spin 
asymmetry B n in the forward kinematics. 
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II. DOUBLY VCS TENSOR BASIS 

In this section, the doubly virtual Compton scatter¬ 
ing (VVCS) tensor will be introduced and the properties 
of the invariant VVCS amplitudes will be investigated. 
We start by reviewing existing developments on invariant 
Compton tensors. A general Lorentz and gauge-invariant 
VVCS tensor was constructed by Tarrach [l£| with the 
corresponding amplitudes free from kinematical singular¬ 
ities and constraints. However, it appears problematic 
to use this tensor basis for practical calculations. There¬ 
fore, we turn to the method used by Prange for RCS 
El and extended to VCS (with only one virtual pho¬ 
ton) by Berg and Lindner [j 2 ■ We consider a scattering 
process 7*(gi) + N(p) —> 7*(g 2 ) + A r (p') and define the 
four-vectors: 


P = 

k = 

L = 


p + p' 

2 

9 i + 92 


91-92 P ~P 


2 2 (1) 

The method of Ref. m for RCS amounts to defining two 
orthogonal vectors, 


P ■ K - 

P'p = pp _ i _il/fM 

k 2 

= E n<*(hp> a Kp L ^, 


( 2 ) 


where we use the convention £ 0123 = — 1 . The four vec¬ 
tors P\n,L'K form an orthogonal basis. Out of these 
vectors, orthogonal and explicitly gauge invariant ten¬ 
sors can be constructed. The six independent Lorentz- 
invariant structures which are needed to describe RCS in 
the most general case are mi: 


Mrcs 


= N‘ 


r pip pw 


(Hi + fCB 2 ) + 


\ P ' 2 

P'»n v - n^P'” 


-(83 + fCBi) 


+ 


We note that 

P^n" -n^P"' 

P'^n-' 


P , 2 n 2 

P’^n" + n^P’ v 
P , 2 n 2 

the terms ~ 


« 75 -B 7 


IB e } N. 


P^n"+n l *P"' 
- pTFpi - *75 


and 

1/1 are ruled out by time-reversal invariance. 
If we allow for one of the photons (normally, the initial 
one) to be virtual, this restriction is relaxed because the 
initial and final states are not identical anymore. The 
method should be modified as follows G 3 : since the con¬ 
dition L ■ K = 0 is no longer valid, we introduce 


L' IJ = 

pip = 

= 


K 2 

p ■ k p t ' 

pv - = — f { P - L 


K 2 

e^P' a kpL' 


L ' 2 


which are orthogonal to each other and to the vector K. 
Furthermore, in order to take into account the longitudi¬ 
nal polarization of the virtual photon, a gauge invariant 
vector is introduced, 


k /fj = k>‘ - 

9 i ■ L> 


( 5 ) 


which obeys qi ■ k' = 0 . With these modifications, the 
most general VCS tensor basis takes the form 0: 


f P'^ P' v 

Kcs = N'l-M-fr+jCB*) 


\ P ' 2 

P'^n v + nPP' u 
P' 2 n 2 

P'> i n v - nPP' v 


-(83 + jCB«) 


P' 2 n 2 
k’l’p’v 


(*75-85 + 1/1B0) 

(*75-87 + i/iBk) 


+ - 


P' 2 K 2 

K'»n v 

i 2 k 2 


(89 + fCB 10 ) 

(*75811 + 1/1B12) 1 V, 


( 6 ) 


The Compton tensor in this form, for the case where 
both photons are virtual does not exist in the literature. 
We follow here the scheme proposed in mj and generalize 
the invariant Compton tensor for the VVCS case. For 
this, we introduce another gauge invariant vector with 
respect to the outgoing virtual photon (72, 


k" v = k" — -k—^-L' u . 

92 • L' 


( 7 ) 


We finally obtain the 18 different structures that con¬ 
tain all the information about Compton scattering with 
two virtual photons: 


Kvcs = N‘ 


( p'P p'v 

P'V + 


-(83 + ^84) 


+ 


(775 i?5 + ifiBe) 


( 3 ) 


P' 2 n 2 
P'»n u - n^P' 1 ' , „ 

+-^2-(^587 + iB 8 ) 

fpP p!v 


+ 

+ 

+ 

+ 

+ 


P l 2 K 2 

k'Pn v 

n 2 K 2 

plpfplv 


P' 2 K 2 
n»k" v 


(89 + fCB 10) 
(*75811 + ^8 12 ) 
(813 + fcB\k) 


n 2 K 2 

k'PK" v 


K 2 


(*75815 + i / iB w ) 

(847 + JCB 1S ) l N. 


( 8 ) 


The invariant amplitudes Bi are functions of the in¬ 
variants: 

Bi = 8j(51,92,91 • 92,8- K). 


( 4 ) 


( 9 ) 
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We next consider the behaviour of the invariant am¬ 
plitudes under photon and nucleon crossing. 


A. Nucleon crossing 


This transformation is a combination of partial parity 


transformation and charge conjugation: "f*(qi) + N(p) 


Y{q2)+N{p') 


l*{qi)+N(-p') -> 7 *(q 2 )+N(-p) 


We note that this is an unphysical transformation since 
only the hadronic part of the kinematics is transformed, 
while the photons remain unchanged. This transforma¬ 
tion, however, relates two reactions which should have 
the same structure. Under nucleon crossing, we have: 


P -P, 
n —> —n, 

k -> k. 

L' -> I/, 

k' -> k\ 
k" -> k", 
k -> c ] kc = -jc, 

•f, —> —ki/c = if,, 
75 ->■ C t 75C = 75- 


( 10 ) 


Hence, nucleon crossing transforms the VVCS tensor 
as 


amplitudes Bi under nucleon crossing: 

Bi{ql,ql,qi ■ q<i,-P ■ K) = +B i {q\,(fc,q 1 - q 2 ,P ■ K) 
fori = 1,3,5,6,7,8,10,14,17 
Bi{q\,ql,qi-q< 2 ,~P ■ K) = -Bi{q\, q\, qi ■ q 2 , P ■ K) 
fori = 2,4,9,11,12,13,15,16,18 

( 12 ) 


B. Photon crossing 


Photon crossing relates two reactions: 


7* ('ll 


N(p) -► 7 *(q 2 ) + N(p') 


7 * (-92) + N(p) 


IK 




Compton amplitude must not be invariant under pho¬ 
ton crossing. It is, however, a relevant symmetry in the 
elastic case q\ = q 2 . Under photon crossing we have: 


k -> 

-k, 

n <-> 

V 

p -► 

Pi 

p' -> 

P\ 

L' -+ 

L\ 

n —> 

-n, 

k' -1 

-k" 

k" -► 

-k'. 


{ pill pm 

___(^ 1 _ ]c n B2 ) 
+^f( N B 3 - k N Bk 

n z 

P ,[l n v + n^P ,l/ N 


(ii) 


pl2 n 2 

P'^rP - rPP' v 


pi2 n 2 


(i 75 n B 5 + pB 6 ) 
(i7s N Br + i N B 8 ) 


J{'l‘ pm 

P' 2 K 2 


( iv p 9 - k N B 10 ) 


K'^rP 

n 2 K 2 

pi tik" 1 ' 

p ,2 k 2 


(i^Bu + i/i N B l2 ) 


{ N B X 3 - k N Bi 4 ) 


n 2 I\ 2 

k'i j k" v 

K 2 


(*75 -B 15 + i 1 Biq) 


( n B 17 - k B 18 ) } N'. 


In the above equation, N Bi stands for the nucleon 
crossing transform of the amplitude Bi. For the invariant 
amplitudes, nucleon crossing amounts to the substitu¬ 
tion P ■ K —> —P ■ K, while the other arguments remain 
unchanged. We next list the behavior of the invariant 


Hence, photon crossing transforms the VVCS tensor 


as 


'Kvcs = & 


, f P'i' P' 1 ' 
\ P’ 2 


CA - k 7 B 2 ) 


(14) 


+ ^-(TBa - k^Bk 


n* 


Vn" (i75 ^ 5 “ * 1Ba) 

P'"rP - n^P'" 


P' 2 n 2 
pi iik" v 


(*75 ^ B 7 - /, J Bs ) 


P' 2 K 2 

rPK" v 

n 2 K 2 

k'li pm 

p ,2 k 2 
k'"rP 


CB 9 - k^Bi 0 ) 
(*75 7 Pn + i/i^Bi-z) 

rBn - k j B 14 ) 


n 2 K 2 

k'> j k" v 


K 2 


(*75 7 Pl5 + l/l^Bi g) 

( 7 P 17 - k 7 B ls )\ N'. 


In the above equation, 1 B i stands for the nucleon cross¬ 
ing transform of the amplitude Bi. Apart from the sub¬ 
stitution P • I\ —> —P • k, photon crossing interchanges 
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the virtualities of the photons as arguments of the invari¬ 


ant amplitudes, q\ <-> 
under photon crossing: 


Bi (q 2 , 

01,01 

• 02, 

-P 

■K) 





for i 

Bi(ql , 

01,01 

’ 02, 

P 

■ k) 





for i 

09(02. 

01,01 

■ 02, 

-P 

■k) 

010 (02. 

01,01 

’ 02, 

-P 

■ k) 

013(02, 

01,01 

’ 02, 

P 

■ k) 

014 (02, 

01,01 

’ 02, 

P 

■ k) 

011 (02, 

01,01 

■ 02, 

P 

■ k) 

012(02, 

01,01 

’ 02, 

-P 

■ k) 

015(02, 

01,01 

’ 02, 

-P 

■ k) 

016(02, 

01,01 

’ 02, 

P 

■ k) 


We find following behaviour 

= +Bi(ql,ql,qi- q 2 ,P ■ K) 
= 1,3,6,7,17 
= ~Bi{ql,q%,qi ■ q 2 ,P ■ K) 
= 2,4,5,8,18 

= -Bu(ql,q%,qi -qi,P- K) 

= B 14 (ql,q%,q 1 -q 2 ,P-K) 

= ~B 9 (qj,ql,qi ■ q 2 ,P ■ K) 
= B 10 (ql,q%, qi ■ q 2 ,P ■ K) 

= B 15 (ql,ql, qi - q 2 ,P ■ K) 

= ~Bie{qt,q 2 ,qi ■ q 2 ,P ■ K) 
= B 11 (ql,ql,q 1 ■ q 2 , P ■ K) 

= ~Bi 2 (ql,q%,qi-q 2 ,P-K) 

(15) 


We now combine photon and nucleon crossing, which 
corresponds to applying a CP-transformation to the full 
reaction. Combining Eqs. GH, we obtain 


-05,8, Bg — B 13, Pro ^ 014, (16) 

011+015,012 + 016 ~ (ql — q%) 2n+1 1 

01,2,3,4,6,7,17,18, Bg + P 13 , B w + P 14 , 

0n 0i5, 0 i 2 0i6 ~ (ql ± q 2 ) 2n > 

with n = 0,1, 2,... 1 

As a consequence, the VVCS tensor simplifies in the 
elastic case (EVVCS), qf = q 2 . 


m EVVCS 


( p'^p'v 

*'{-p 2 - (Bl +^B 2 ) 

+ ^-( 0 3 + jCBi) 
n z 

P'^tx v + n»P ,v 


P r2 n 2 

P , ' 1 n v - n»P' v 

-rvT - 2 -*7507 

P' 2 n z 

//'/' ptv _|_ pl^f'IIV 

p r2 k 2 

k' ti n v - n»K" v 


(Bg + jCB 10 ) 


n 2 I< 2 

K'»n v + n»K" v 


*75011 


+ - 


n 2 K 2 
K'i J k" v 


->12 


K 2 


(0i 7 + fCB i 8 ) } N. (17) 


1 We do not expect any negative powers of q\ — q 2 since it would 
lead to an unphysical singularity. 


We finally investigate the RCS limit, i.e. let the virtu¬ 
alities of the photons go to zero. In this case, 

K" v <£. (18) 

In an observable, the Compton tensor is contracted 
with the photon polarization vectors, £\(qi) ■ K' —> 
£\(qi) ■ qi = 0. Though the amplitudes 0g, 10 , 11 , 12 , 17 ,is 
should not necessarily vanish in this limit, the corre¬ 
sponding terms do not contribute to any observables. 
Therefore, at the real photon point we recover the old 
RCS tensor of Prange of Eq. 0 . 

It should be stressed that doubly VCS tensors for in¬ 
elastic (ql 7 ^ ql) and elastic (ql = ql) cases in the form 
of Eqs. 181171) do not exist in the literature. 


III. APPLICATION: BEAM NORMAL SPIN 
ASYMMETRY 


We consider elastic electron-proton scattering e~(k) + 
N(p) —> e~(k') + N(p’), with the electron beam polarized 
normally to the scattering plane with the polarization 
vector defined as 





k x k' 

| k x k'\ 


(19) 


The single spin asymmetry defined as B n = 
is refered to as beam normal spin asymmetry. It was 
shown that such an asymmetry, to leading order in the 
electromagnetic coupling constant a em ~ 1/137, is given 
by 


Espins 2Im(A4C 7 AbsA4 27 ) 
Espins IX17I 2 


The absorptive part of the two-photon exchange (TPE) 
graph is given by 


AbsAd 27 = e 


|fci| 2 d|fci|dn fcl _, 


20 i ( 2 tt ) 3 


u’^ v (h + 771 ) 7 ^ 7 * 


w Abs M^(w 2 ,t,Qi,Q 2 2 ), (21) 


QIQ 


where is the VVCS tensor. Furthermore, the asym¬ 
metry can be written as 


0 „ 


e 2 t 

D(s, t) 


|fc 1 | 2 d|£ 1 |df 4 1 1 

2E 1 (2nf Q 2 Q 2 


lmL a ^H a ^(22) 


The leptonic tensor L alll , is given by 
L atlv = Tr(/f' + 771 ) 77,(^1 + m) 7 M 7 5 ^(^ + m) 7 a , (23) 


while the hadronic tensor is defined as 


H AbsM /il/ (A r 'r a A r )*, (24) 

spins 
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with T“ = Gm{~ t)"f a — F 2 {—■ In the following, we 
will use the hadronic tensor in the form introduced in the 
previous section for the most general inelastic case. The 
structure of the tensor of Eq. © can be represented as 

■ N'N + ■ N'fCN 

+ C^ ■ N'i^N + D^ u ■ N'lfN. (25) 

The tensors (A, B, C, D)^ v are defined in Appendix. 
The sum over spins in the hadronic tensor can be readily 
performed with the result: 


H= 8MP 1 

+ 8 L 2 G 


PK 

G e A^ + —Ei 


LK 


(26) 


M 


I< a - —L a B^ + n a D^ u 


Notice that the tensor C appears with 75 , and does 
not contribute to this observable. We list here the re¬ 
sult of the tensor contraction and refer the reader to Ap¬ 
pendix for the details of the calculation: 


La^H 01 ^ = 64 imM(n-£) 

( f T2n2m2 

X 


S Ge (b 1 + b 3 -F1 M s A 
I V K 2 (L 12 ) 2 J 

PK L 2 Q\Ql \ 

AP^PK) f g* ) 

P' 2 K 2 \ M 


(27) 


Q! 

0 

n z 


PK - pkiJ m, £] { a * B ° +s w F ' B ™ 


Ql 

n‘ 


t\ \ pk+ pk^ l2 ] (r r + pk FR 

? [ ph + ph ~k^\ l ^ 13 + ~M FlB 

8 imL 2 G M ^l 

B S - (P'K)B 2 + 9l±M iQ 2 Bl0 + q 2 b 14 )I 
8 

n z 


\lk)b s+ mA (Q 

V ' 1 S' 2 r /2 ^ 


-b 4 


- l 2 g 


■mB 6 


For forward kinematics, the virtualities of the photons 
are very close to each other such that the contributions of 
the terms ~ {Q\ — Q 2 ) n are suppressed by powers of t/s. 
In this case, we can use realtions between the amplitudes 
Bi established in the previous sections. Denoting Q\ = 


Q\ = Q 2 , we obtain: 

L aiJ ,vHzv VCS = 64 imM(n-£) 


(28) 


x "I Ge \ B\ + B 3 — 


Q 4 


I\ 2 L’ 2 


Bn 


Fi \ B 2 + B 4 — 


PK 

~M 

2(P'k 1 )(PK) 


Q 4 


K 2 L' 2 


-B 


18 


P’ 2 K 2 


' PK 

G e B\ A -j-j-F\B 2 


+ 


2 Q 2 


PK + PK- 


L 


2 1 


K 2 


G e Bq 


PK 


■F 1 B 


10 


-128imL 2 G„(n{) J 

-128i m KAAAAl 


x | M 3 G e 
- G m 


, t PK „ 

3+ AT 54 


L 2 B 6 — (PK)S—Bi 2 


A further simplification can be obtained in the quasi- 
RCS approximation. It amounts to considering a soft in¬ 
termediate electron. This kinematical region is enhanced 
by large logarithms ~ ln(— t/m 2 ). This approximation 
corresponds to taking k\ = 0 and Q\ = Q 2 = 0 in the 
numerator. 


= 64 imM(n-£) 


(29) 


. PK 

x \ Ge (Bi + B 3 ) + (B 2 + B 4 ) 


IV. CALCULATION OF B n IN THE FORWARD 
KINEMATICS 


In this section, we are going to apply the formalism 
developed in the previous section to a calculation with 
Regge kinematics. Recent calculations have shown that 
in the forward kinematics, a logarithmic enhancement 
takes place [l3j, but the contribution of the In 2 (—t/m 2 ) 
term comes with photon helicity flip am plit ude of RCS 
and its contribution is negligibly small Il4 0 . Van¬ 
ishing of the RCS contribution at the exact quasi-RCS 
point requires that we take into account the full k\ and 
Q\ 2 dependence of the hadronic tensor. The only sim¬ 
plification which arises due to the forward kinematics is 
based on the observation that integrals containing pow¬ 
ers of (Ql — Q 2 ) are suppressed by powers of the (small) 
momentum transfer t. It can be shown, in fact, that 


\k 4 \ 2 d\k 4 \<m kl Q\-Q\ 
2Ei(27t) 3 Q\Q\ 


f \ki\ 2 d\ki\dCl kl (Q\ - Ql) 2 
J 2E 4 (2tt) 3 Q\Q\ 


0(t/s), 


(30) 
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and we refer the reader to the Appendix for details of 
the calculation. Therefore, in the Regge regime, elastic 
VVCS should give the leading t contribution, and we can 
use the hadronic tensor in the form of Eq. ra whose con¬ 
tribution to B n is given by Eq. (GEJl- In order to calculate 
the asymmetry B n in these kinematics, we are going to 
relate the invariant amplitudes P; to the helicity ampli¬ 
tudes of doubly VCS. The imaginary parts of these latter 
are related to the total photoabsorption cross section by 
means of the optical theorem. 


where we explicitly took out the electric charge e. We 
next define shorthands, 

Tll = IoA'jOA 
T T t = T± A';±A 
TtL = T± A' ; oA 

Tlt = Toa' ; ±a- (36) 

With these definitions, we obtain: 


A. Helicity amplitudes of VVCS with Q\ = Q\ 


For this calculation, we should specify the reference 
frame in which we work, since the helicity amplitudes 
are frame dependent. This frame dependence is then 
eliminated with use of the invariant amplitudes. For the 
calculation, the Breit frame will be used 2 . It is defined 
by p + p' = 0. Furthermore, 


0), (31) 

= (w,0,0 ,q), 

V* = ( 0 ,y, 0 , 0 ), 

91 = (^^, 0 ,g), 

92 = (w,-y,0,g), 

where A = > 0. We next define the polarization 

vectors of the virtual photons with the photon helicity 
Ai = ±1,0 referring to the initial photon and A 2 = ±1,0 
to the final photon: 

e?(Ai =± 1 ) 
e?(Ai = 0 ) 

e£(A 2 = ±l) 
e 2 (A 2 = 0 ) 

(34) 


^( 0 , «i’ iAl ’ 291 )’ 

= v^(‘ , W‘ a )’ (32) 

= “ 7 i( 0 ’P iA2 '%)' 

= v^(‘ , W‘ s )’ (33) 


where <ft denotes the magnitude of the photon three- 
vector, qi = | g! | = \q%\ = \J q 2 — t/ 4. The Compton 
helicity amplitudes are defined as 

TxzX'-mx = e 2 £^(Ai)(4(A 2 ))* • M m „, (35) 


2 Note that the 7 *p c.m. frame cannot be used here, since in 
this frame it is impossible to fix the electrons’ kinematics to the 
experimental one 


jT ll = 


4 Q 2 


t + 4 Q 




r<z2 (Pi+^P 2 ) - ^(Pit + A'Pis) 


9i 


9i 


2lo 


V^ 2 - y 2 


(B 9 + ]KB W ) 


Nx, (37) 


-T T t = ~N' X , 


t 


(Pr+^P 2 ) 


[qft + 4 Q 2 
-ArA 2 (B 3 + f(B 4 ) 

4Q 2 2 uj 

+ t + 4Q 2 v / M 2_| g2 


(Bg + fxB w ) 


16 Q 4 g 2 

t(t + 4 Q 2 ) qf 

— (Ai + A 2 ) 


(Bn + fCBig) 


yjM* - 


-Ik 


16 Q 2 


lit 


+ (A 2 — Ai) 


sJm 2 - i9i 

2 a; „ 16Q 2 


t(t + 4 Q 2 ) 


75 


B 


12 


-Ik 


y/ M 2 - \q\ + 4| 5 2 ) 


B, 


(38) 


Nx, 


jTtl = — 


(39) 

4W -‘“ 2 1 -(B, + /CB 10 ) 


<7 2 (t + 4Q 2 ) _ 

4Q 2 ~ 2 


t + 4Q 2 qf 


(Bn + fiB ig) 


N x , 


vWa 2 

(AP - i)A q Ql 


NyfaBr - if,Bk)Nx 


4y / 2Q^q;A2A^(75Pii - ijB 12 )Nx 
\J M 2 — jAqqi(t + 4Q 2 ) 
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VWuq R , \ t 

1 LT A q\ V [t + 4Q 2 


f _(B! + fCB 2 ) (40) 

4QV - fw 2 1 


q 2 (t + 4Q 2 ) ^/ M 2 _ | 
4Q 2 g 


-(B 9 + /STB 10 ) 


i + 4Q 2 g 2 


(Sl 7 + fiB is) 


Aa, 


vWAi 


(M 2 - l)Aq qi 


N' x , ( 75 -B 7 + i),B 6 )N\ 


4 v / 2Q 2 wAi 


M 2 - {A qqi (t + 4Q 2 ) 


N'^^Bn + ii/iB 12 )N\, 


We next note that the amplitudes B 7 and Bn do not 
contribute to the beam asymmetry. As can be seen 
from Eq. ®, the amplitudes Bq and B 12 depend lin¬ 
early on the photon helicities and therefore drop out of 
the ’’cross section combination” of the helicity ampli¬ 
tudes, (Tia'jIA + T_ ia' ; -ia)- Instead, they obtain their 
contribution from the GDH sum rule-like combination 
(Tiv; U — T-ia';-ia) which vanishes in the considered 
kinematical regime according to Pomeranchuk’s theorem. 
In previous work, the contribution of the photon helicity- 
flip Compton amplitudes was studied for the beam asym¬ 
metry [ 13 , [l|j ■ In spite of double-logarithmic enhance¬ 
ment due to the quasi-RCS kinematics, their contribution 
to the asymmetry B n was found to be negligibly small 
for forward kinematics, and will be neglected in the fol¬ 
lowing. 

The longitudinal-transverse amplitudes are in general 
non-zero. However, it is only the combination Ttl—Tlt , 
with the conserved photon helicity Ai = A 2 that con¬ 
tributes to the amplitudes 131 , 2 , 3 , 4 , 9 , 10 , 17,18 relevant for 
the asymmetry. As a microscopic calculation shows (see, 
for instance, |TJ) this difference vanishes for diffractive 
kinematics. Thus, the only non-zero helicity combina¬ 
tions are E T t = ^(7ia' ; ia + T_ia';_ia) and T LL . Their 
imaginary parts are related to the transverse and longi¬ 
tudinal cross sections, respectively, which contain all the 
physical information on forward doubly virtual Compton 
scattering with an unpolarized target. For the remaining 
amplitudes Bi , we define shorthands, for example Bi for 
N' x ,(Bi + ]fCB 2 )N\, and find the expressions as functions 
of the helicity amplitudes: 


Working out the spinors, we get in the Breit frame: 


N'N = 2y AI 2 — ^2A5\-\>, 

N'fCN = 2Mui2\5 X -\> + q A2\S X y. (42) 


We see that in the forward limit t = 0, only nucleon 
helicity-flip amplitudes survive.This should not mislead 
the reader since the reason for this helicity flipping is the 
use of Breit frame kinematics for the nucleons: the initial 
and final nucleon three-momenta are opposite by defini¬ 
tion, therefore unchanged spin projection onto the x-axis 
will correspond to opposite helicities. Furthermore, the 
optical theorem relates the imaginary parts of the for¬ 
ward amplitudes to the corresponding cross sections, 

Im^T^i.ii +T_ 1 _i._ 1 i) = 2yJ M 2 — t/Auiar 

ImT 0 _ 1 ;0 i = 2\JM 2 - t/4ua L (A3) 


If the Callan-Gross relation holds, the longitudi¬ 
nal cross section is related to the transvers one as 
Q 2 

<tl = ^• Inserting these expressions into the ten¬ 


sor contraction of Eq. |28|l and noting that Abs M^ v = 
— 2ImM / “', we obtain: 


L atlv H^ v vcs = = -128 imMGMn • 0 


2uE 1 Q 2 Q 2 Q 4 


q\ 


4 , 

Qi i e 


q q 
128i »«VC« |,,tl) f' fl 4 




(44) 


It can now be seen that spurious singularities ~ t+iQ 2 
appearing in the individual amplitudes in Eq. m cancel 
in the observable, as they should. 

We next perform the integration over the electron 
phase space. We remind the reader that in the choosen 
kinematics 


n? = ( 0 , 0 , —\Jm 2 — ^ 9 ^, 0 ), 

e = (0,0,1,0), (45) 


e 2 Bi 

e 2 B fl 


UJ 


t 


3 

e 2 B 9 

r 24Q2 B 

e — B ir 


q\ t + 4Q 2 
—E tt , 

~ (7 2 

(PK)« 


E tt + 


4Q 2 


ql t + 4Q 2 


Tll, 


t 


qft + AQ 
2 4 Q 2 


r(E tt — Tll), (41) 


q{ t + AQ‘ 


LjTT — 


t 


q{t + 4 Q 


;T L l- 


therefore 


|fci| 2 d|fci|dH fel 1 

2Ei(2tt ) 3 Q\Ql 


= -128mMG E \lM 2 - ~^<tt 


2 Eq - 2J 2 - 2/2 + h ~ 


2 M 2 


M 2 — | 


J 25 


(46) 































where we have used the notation: 


the final result for normal beam spin asymmetry: 


Jf 

h 

h 


j r 


|fci| 2 d|fci|dn fcl fcf 


2 J B 1 (2^) 3 Q 2 Q 2 


|fci| 2 d|fci|df2 fcl 1 

2£'i(27r) 3 2 

|^| 2 d|^|dn fcl (Q? + QD 


1 1 

Ql + Ql 

2 
1 


2 Bi (2^) 3 AQlQl (k - k,f 


|fci| 2 d|fei kiki 
2E ± (2tt) 3 Q 2 Q 2 


(47) 


-Bn = — 


m 




G e 


4tt 2 (1 + t) 2 tG 2 m + eG 2 E 


. , t E 2 

Inf-) -l + _ln_ 


m J 
t 


AE 2 


777/ 

M 2 


ln W + 1+ 2 s 


( 52 ) 


V. RESULTS 


Furthermore, the coefficient J 25 originates from decom¬ 
posing the tensor integral according to na, 

J' 1 " = J 21 P»P V + J 22 K ll K 1 ' + J 2S {P ll K 1 ' + K^P”) 

+ + J 25 ^- (48) 


The integrals I 3 and I 2 entering Eq. ll47ll were calcu¬ 
lated in the previous work 0 . We obtain for the first 
three integrals 




_ B f Ethr \/ t Ethr 


-AttH ( E 


E 


(49) 


+ 1 


- l n (1 _ 


E J 


1-2 


■ 


B = 


1 Ethr 2 Ethr 

In- 


4tt 2 B 

1 

47T 2 
1 

87T 2 
1 

87 r 2 


m 


E 

1 + %) 

Ethr \ 


1 - 


B J 


b y 


In 1 + 


ln 1 — 


Ethr \ 

y 

Ethr \ 
Ethr \ 

J 


(50) 


and for the tensor integrals: 


too 

J 2 


J25 — 


E 2 thr 

— 8 n 2 t 

E 2 

— 87 r 2 t 
1 


ln 


QE, 


thr 


E 


1 - 


Elhr 


E. 


thr 


B 2 


In 1 - 


—) 

Ethr J 

Ethr 


E 


32t r 2 B 2 + t /4 
B 2 


i„/5+ 


1 


2 B 


B 

Ethr 2 


(51) 


+ 


E 2 thr 


- 1 In 


Ethr \ 


B, 


thr 


327 x 2 [pk) 


In the above formulas, E t h r denotes the upper limit in 
the integral over the electron energy, which corresponds 

to threshold production of the pion, E t hr = 

in the Breit frame. We refer the reader to the Appendix 
for details of the calculation. We can finally write down 


The model used in the previous section for the calcu¬ 
lation of the B n is based on the optical theorem, Callan- 
Cross relation between the longitudinal and transverse 
partial cross sections, and the neglection of the mismatch 
between the virtualities of the virtual photons inside the 
loop. This latter assumption restricts the applicability of 
the above formula to low values of the momentum trans¬ 
fer. To estimate the error introduced by neglecting the 

terms ~ fo/A under the integral over the electron 
phase space, we note that the integrals entering the final 
result of Eq. (I521l have the order ~ f _1 (I 3 , J 2 °) and t° 

_ q! \ 2 

(I 2 , T}, J 2 h )■ As it was noticed before, terms ~ 

contribute to the order t 1 . Therefore, they contribute at 
relative order of t 2 . Thus, for low values of t , our full 
result of Eq. (15^1) is consistent with the approximation 
used. 

We next study the leading f-dependence of the beam 
normal spin asymmetry. It is given by 


B„ 


niy/—t<7T Ge 

471-2 tG m + £ G 2 e 



1 (53) 


This result is in exact agreement with the result of 0 . 
The approach of 0 consisted of taking only one ampli¬ 
tude for elastic doubly VCS, the one which survives in 
the exact forward RCS. It therefore amounts to neglect¬ 
ing not only the additional t dependence of the Compton 
amplitude, but also the dependence on the photon vir¬ 
tualities. The present calculation, on the contrary, uses 
the complete set of invariant amplitudes for the case of 
equal virtualities of the incoming and outgoing photons, 
and does not apriori set any kinematical variable to zero. 
The reason for this agreement is that every factor ~ Q\ 2 
in the numerator actually leads to an extra suppression 
in t, as can be seen, if comparing the integrals I 2 to Iq 
and B(* to If, respectively. The integrals are given in the 
Appendix. 

If using as input the exact forward value of the imag¬ 
inary part of the Compton amplitude, to go to finite 
values of t, one needs a phenomenological information 
in order to consistently describe the off-forward Comp¬ 
ton scattering. We use the same approach as 0 . It is 
based on the phenomenological exponential fit of the in¬ 
dependence of the Compton differential cross section for 
















































values of t from 0 through -0.8 GeV 2 . Since ^ ~ 
one gets 
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crr(t) = ctt( 0)e Bt ^ 2 , (54) 

with B ss 5.2 — 8.6 GeV -2 . For the energy dependence 
of er T , we use the phenomenological fit of Ref. djj- We 
present our results in Fig. UJ As seen from the figure, 



FIG. 1: The results for B n as function of the momentum 
transfer Q 2 are shown for four values of the lab beam energy: 
3 GeV (solid lines), 6 GeV (dashed lines), 12 GeV (dotted 
lines), and 45 GeV (dash-dotted lines). The thick lines corre¬ 
spond to the full f-dependence of Eq. (1.521 . while the thin lines 
correspond to the leading t behaviour of Eq. 1531 . 

the more accurate account on the t dependence leads to 
faster decreasing of the asymmetry with increasing Q 2 . 
We predict therefore lower negative values for B n than 
the authors of [TdJ. 

In Fig. [3 results for B n on the neutron target are 
shown. The asymmetry is only about twice smaller for 
the neutron than for the proton, which is due to the 
overall factor of . Since also the cross section in 

T ^M+ ELr E 

the denominator is much smaller for the neutron target at 
low Q 2 than for the proton, it means that a measurement 
of the cross section difference in the numerator of the B n 
would be a very challenging task for the neutron target. 


VI. SUMMARY 

In summary, we constructed the invariant basis for 
doubly virtual Compton scattering (VVCS). In the gen¬ 
eral case, VVCS is described by means of 18 independent 



Q 2 (GeV 2 ) 


FIG. 2: The results for B n as function of the momentum 
transfer Q 2 for the neutron target. Notation as in Fig. 0 


invariant amplitudes which are functions of the invari¬ 
ants PK, t , fZi, <72- The proposed form of the VVCS ten¬ 
sor does not exist in the literature. This form is based on 
Lorentz and gauge invariance and consists of constructing 
9 orthogonal basis tensors, which are useful for practical 
calculations. We furthermore investigated the behaviour 
of the VVCS amplitudes under photon and nucleon cross¬ 
ing. Basing on this crossing behaviour, it is possible to 
obtain the correct form of the VVCS basis in the case 
of elastic VVCS (i.e., qf = q 2 ) where only 12 ampli¬ 
tudes survive. The special case of elastic VVCS is real¬ 
ized in the case of the calculation of the beam normal 
spin asymmetry B n with forward kinematics. We pro¬ 
vide the calculation for this observable with these kine¬ 
matics. Since in the exact forward limit this observable 
vanishes, one has to use the full amplitude for this calcu¬ 
lation. We project the Compton helicity amplitudes onto 
the invariant basis and use as input the optical theorem. 
It relates the imaginary part of the helicity amplitudes 
+T , _ii _ 1 i) and T 0 i Q i to the transverse a t 
and the longitudinal a L cross sections, respectively. The 
longitudinal cross section is then related to the transverse 
one by means of the Callan-Cross relation. The resulting 
asymmetry is of the order of ss —5 ppm and is practically 
independent of the electron energy. 

The presented formalism may be furthermore extended 
to the two photon exchange contributions to other ob¬ 
servables. This work was supported by the US Depart¬ 
ment of Energy contract DOE-FG02-05ER 41361. 
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VII. APPENDIX 


The leptonic tensor 


A. Projection technique 

In a real calculation, one needs to calculate the single 
amplitudes contributions from, for instance, a Feynman 
diagram. For this, one can make use of the fact that all 
the basis tensors are orthogonal. For example, 


pf/j, pfiy 


B 1 + fCTh 


N, 


(55) 


and so on. We find following useful relations: 

„ , PK L 2 

p'^ _ p* _ 


L ' 2 = L 2 - 


K 2 L ' 2 
(LK ) 2 


K 2 

n 2 = ~P' 2 K 2 L ' 2 
K 2 

K ' 2 = J h Q2 

K 2 

k " 2 = -%Ql 


K 2 

U 2 


k' ■ k" = ^ (l 2 - k 2 ^j , 


(56) 


and furthermore, 

P' ■ K = P' ■ L' = P' ■ L = P' ■ Ik = P' ■ K" = P' -n 
= n-k = n-L' = n- L = n-k' = n- K" = 0; 
qi ■ L' = -q 2 ■ L' = L ■ L' = L ' 2 . (57) 


B 


Lain, = Tr(^' + m) 7 „(^i + m) 7 ^ 75 ^ + m) -y a (62) 


can be conveniently decomposed into two parts which are 
symmetric and antisymmetric in the pair of indices (fM ^): 


Leptonic and hadronic tensors 

adronic tensors quoted in Eq. < 1251 ) are given by 
piii pm k't‘P w 

p ' 2 Bl + ~k~ B:i + p' 2 k 2 Bl> 

piII p/iV pill pirn 

+ - B 13 + -=- B 17 (58) 

p/2 J p 2 v ' 

piii pm n fJ, n 1 ' k'»P’ v 

' = ~r ^ B2 + — Bi 

p'nk" v „ k'»k" v „ 

+ ----^14 ~ -^18 (59) 

p'2j{2 J(2 10 v ’ 

P'V + n^P'" D P^n v - n^P ,v D 
Pikf- 35 + P^k , 2 Bl 

k'^n v „ K' w n p „ 

-F - ~—Bn -=—S15 ( 60 ) 

n 2 K 2 n 2 K 2 V ’ 

.. + n p P' u P'^nP - n p P' v 

> 12 n 2 P' 2 n 2 B& 

K" v rP 

** + kpW Ble (61) 


a™ = -8 im{g liU e K p Ta L K k p t T (63) 

+ {kl) l i£ KV TaL l 'k + (fci)„£ KMTa L K £ T , } 


rASYM 

J CX.^LV 


= - 8 im{£ u £ K , l _ lTa L K K T - ^£ k „ toi L k K t 

T Lct£vilTpfi L £vprck, 


- K a £uiirpC KP + £, a £viiTpklL p 


} ■ (64) 


C. Tensor contractions 


L 2 K a - ( LK)L a B pv ■ L a/n/ = -8 imL 2 {n^) 


2(P'ki) Q 2 + Q\ r n 2 D . ^2 ] 


p /2 

16im(L 2 ) 


■B 2 + 


2 s 2 (nki)( K ^) 


n 


4 ?r 2 

2 (PK)B 4 . 


[QlBio + Q2B14] 


(65) 


n a D^ ■ L atlv = 16 

n z 


x 


( 66 ) 

-L 2 Bq + ( LK)Bs + ^ j \Q\B 12 + Q\Biq\ \ . 

K Z L' Z ' 1 
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P a A^ ■ L aflv = 16 im 


(nki)m) p2L 2 B 


to the integral 

' |fci| 2 rf|fci|df2 fel ( Ql-Ql ) 2 


h = 


2E 1 {2ir) 3 Q\Ql 


= 16 L„L V J^. (69) 




n* 


The first integral was already calculated BUI- We 
next turn to the vector integral 


+ 


+ 


[k'K") 2/v 2 (q 1 L)(K"k 1 ) - (q 2 L)(K'k 1 ) 


A 2 A ' 2 

(■ Pk)L 2 ( (A'fci) 
n 2 ^ K 2 J 

(PK)( qi K) + (P'h)( qi L) 


K 2 L’ 2 


- 1 


B 


10 


(PK)L 2 / (K"k r) 

n 2 1 A '2 


- 1 


(PK)(q 2 K) - (P'h)(q 2 L) 


B\a 




+ 8 *m(n£) |Ai + P 3 + 

L 2 Q 2 Q| 


2(P'k 1 )(PK) 


n* 


Bi 


K 2 (L’ 2 y 


B 17 


, Of 

' o 
n z 

, Q| 
' 2 
n z 


(PK) - ( PK) 


(g 2 A)A 2 


(PA) + (PA) 


K 2 L' 2 _ 
( qi L)L* 


I\ 2 L' 2 


IP 


B 13 


D. Integrals over the electron phase space 


The integrals appearing in the final result are: 


In = 


/f = 


h = 
h = 


J? = 


lAhpdlfcrldn*, 1 

2Ai(2^)3 Q 2 Q 2 

|fci| 2 d|fci|dfl fcl fcf 

2Ai(2^)3 Q 2 Q 2 

|fci| 2 d|fci|dfl fcl 1 r 1 1 

2A 1 (27t) 3 2 [q 2 Q 2 J 

Ik^dlk^d^ {Ql + Qlf 1 

2Ai(27t) 3 4Q 2 Q 2 (k-ki) 2 

\ki\ 2 d\ki\dfl kl k^kl 
2A 1 (27t) 3 Q\Q\ 


Bn 


i r = 


P^thr 


k\dk\ 


dflk 


K 


0 2A 1 (2tt)3 J Kl (fc_fc 1 ) 2( fc /_ fcl )2 

= / 1 P P^ + Iir-K 1 ' . (70) 

It cannot depend on due to the symmetry of Ji 
under interchanging k and k'. To determine these two 
coefficients we have a system of equations, 


k°i 1k + P°hp = A u = 


i 


d 3 ki 

16tt 3 J Q\Q 2 


(71) 


II ik + 4P Klip = 4A / ,/f 


d 3 fci 


= I 2 - 


(67) 


J (2tt) 3 E 1 Q 2 2 ’ 

Using the same approach as for 7 0 (see the Appendix of 
B) , we obtain for J®: 


A u = 


A 


zdz 


— 87 r 2 t 


z 2 - ^ - ^{1 - z) 2 


(72) 


■In 


-2 _ _ 4m 2 ^ _ ^2 | / ~2 m 2 


.2 _ _ 4m 2 ^ _ £^ 2 _ 


2 _ m 2 
4 TJ? 


E f Eth r \/ t Eth r 
—A-kH \ E m E 


E, 


thr 




(73) 

(74) 


Finally, we consider the integral I 2 where we perform 
the angular integration, 


h = 


1 


d 3 k\ 

( 2 tF J E^Ql 

,J thr 


(75) 


1 

8 ^ 


dz In ■ 


1 




After integrating by parts and changing variables, z = 
yr cosh y and y = In t, we arrive at 


( 68 ) 


r 1 A thr 1 2 Ethr 

-'2 — - 7-0 —FT - - 

47 t 2 E m 

1 ( Bthr \ j A thr\ 


(76) 


47T 2 


1 - 


a y 


A 


T 


/q 2 _q 2\2 

Furthermore, the terms ~ , which we have 

Wi W 2 

consistently neglected throughout the calculation, lead 


In these last two integrals, it is important to keep Aj^r. 
unequal to 1 until the end to ensure the convergence of 
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the integral. Solving the system of linear equations for 
the coefficients, we obtain: 


hp = 
I\k = 


1 s-M 2 E thr 2 E 
8 t r 2 M 4 — su E U Q 


1 



APR 

~Q 2 ~ 


hp- 


(77) 


We next consider the integral I 3 and rewrite it as 


h 


2E 1 (2n) 3 


(Ql-Ql) 2 ' 

4 Q\Ql 


1 

(k-k 1 ) 2 ' 

(78) 


We note that the second term in the square brackets 
~ (Q 2 — Q 2 ) 2 is of the order t/s. For the moment, we 
keep this subleading term, and calculate the leading in t 
contribution: 


h — 


\ki\ 2 d\k 4 \dCl kl 


1 


1 

87T 2 

1 

87T 2 

1 

&K 2 


2E 1 (2tt) 3 (k-k!) 2 

f kM , . 1 + z 
/ dz m 


fO 


l-z' ° {t/s) 

1+ ¥) b ( 1+ % : ) 

1 _ In 


0(t/s).( 79) 


As we see, the leading term contributes to B n to 
relative order t 1 . This means that the second term 
should be omitted to be consistent with the approxima¬ 
tion adopted. 

The calculation of the tensor integral is more involved. 
We decompose this tensor into tensorial structures con¬ 
structed out of external momenta 0 , 


2 k! v J% v = 2PK J 24 P^ + 2AT 2 J 22 K^ 
+ 2 [PKK' J + K 2 P>*] J 23 
- 2L 2 J 24 L^ + 2J 25 (K - L)^ 


f \k!\ 2 d\k!\dn kl k j* 

J 2E!(2tt) 3 Q\ 

Rewriting the vector integral B\ as 

\ki\ 2 d\ki\dn kl k J* 


= Bf. 


(83) 


Bf = 


= Bn/ + B 12 k» 


j 2Ei(2n ) 3 Q\ 

= B n P» + B 12 K» + (B 12 - B n )A/ (84) 

we find following relations between J 2 i and Bp. 

2 PKJ 21 + 2K 2 J 2 3 = B n 
2K 2 J 22 + 2PK J 23 + 2 J 25 = B 12 

2L 2 J 24 + 2 J 2 5 = Bn — Bi 2 . (85) 


To proceed, we calculate the appearing vector integral. 
The coefficients can be found as following: 


2 fc M Bf = 2pkB u = C 0 , ( 86 ) 

while the 0-component B 3 = p°Bn + k 0 Bi- 2 can be cal¬ 
culated directly: 


B° 


\ki\ 2 d\ki\dn kl 1 


167T 3 
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167T 2 


El r 

167r 2 B 


Q? 

2 + 


In ■ 


In 


2 — 

l 

m 2 


2Ethr 


(87) 


_ f \ki\ 2 d\ki\dCl kl Kkl 
2 J 2Ei(2tt) 3 Q 2 Ql 

= JnP [l P v + hiK^K" + J 23 (P /i A"' + P V K») 

+ J M I/ t L v + ■ (80) 

The final result contains In order to calculate the coef¬ 
ficients J 2 j , we contract the tensor integral with and 
external momenta: 


J 2 


A 4 


P 2 J 21 + A' 2 J 22 + 2 PKJ-23 + A 2 J 24 + 4J 25 


2 f \ki\ 2 d\ki\dfl kl 1 

J 2Ei(2tt) 3 Q\Ql 


(81) 


2k,2k' v J^ 


= (2PK) 2 J 21 + 4(AT 2 ) 2 J 22 
+ 8K 2 PKJ 23 -A(L 2 ) 2 J 24 + 8K 2 J 25 


\k!\ 2 d\k!\dn kl 

2Bi(2tt) 3 


= C 0 = 


87T 2 


(82) 


The fifth independant equation for five unknown co¬ 
efficients J 2 j is obtained by the calculation of the J 30 
component of the tensor integral. The calculation fol¬ 
lows the same steps as that for I 3 given in the Appendix 
of M- We obtain: 


rOO 

J 2 


\k!\ 2 d\k 4 \dn kl 1 

-Ei 


E 2 


167T 3 

k+h.T 


Q\Ql 

z 2 dz 


-8ir 2 t 


z* - 


‘-(l-z) 2 


In 


z + 


- ^f(l-zf 
- 4=4(1 - z) 2 - z 


( 88 ) 


Due to two powers of 2 = k 4 /k in the numerator, the 
contribution of small values of the intermediate electron 
momenta is negligible, and the square root can be Taylor- 
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expanded since m 2 /\t\ « 1. This leads to 


700 _ 
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E 2 thr 

— 87 i 2 t. 
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(89) 


Solving the system of linear equations for the coeffi¬ 
cients J 2 j , we find 


J‘25 — 


p2 

p/2 
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too _ E 
2 2 K 2 

E\ hr 


32tt 2 E 2 
E 2 

+ ' 
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7 4(pfc) 
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(90) 


327T 2 (pfc) 


In the latter equation, the vector P is defined as 
p/M = pi 1 — with P' 2 = P 2 - where 

K* = \{k + wy. 


Finally, the integral / 4 is given by 


h 


2\2 


f Ifcrpdlfcrldfifc, {Ql-QD 

J 2£i(2tt) 3 4 5 6 7 8 Q 2 Q 2 

8L 2 (2L 2 J 24 + 2 J 25 ) = 2t{B u - 
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(91) 
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